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1. Introduction 

One-photon double ionization of He has been a subject of intense theoretical and 
experimental studies over the past decade. As a result, a complete understanding of 
this reaction has emerged both in terms of underlying physical processes and the precise 
knowledge of magnitudes and shapes of the cross-sections over a wide range of photon 
energies (Briggs k Schmidt 2000, King k Avaldi 2000, Avaldi k Huetz 2005). 

Two-photon double ionization (TPDI) of He represents a new level of complexity 
and brings new challenges, both to theory and experiment. Extremely intense VUV 
radiation is needed to study TPDI of He. Only very recently, such a radiation has become 
available from the free-electron laser sources (Laarmann et al 2005) and high harmonics 
generation (Nabekawa et al 2005). Dynamics of the TPDI is richer as compared to 
one-photon double ionization because of interplay of the S and D continua. Also, the 
target description becomes more involving since the knowledge of the intermediate state 
as well as the initial and the final states is needed. More importantly, the target states 
can be modified by the strong laser field. 

There have been several reported calculations of the total TPDI cross-section of He 
at various photon energies (Pindzola k Robicheaux 1998, Nikolopoulos k Lambropoulos 
2001, Parker et al 2001, Mercouris et al 2001, Colgan k Pindzola 2002, Feng k van der 
Hart 2003, Piraux et al 2003, Hu et al 2005). Although numerical values of the cross- 
sections vary depending on the theoretical model and assumed characteristics of the 
laser field, there is some consensus between several calculations (see Hu et al (2005) 
for detail). It has also been well established that a non-perturbative treatment of the 
electron-photon interaction is needed to obtain accurate total cross-sections. 

As compared to the total cross- sect ion, the present knowledge of the differential 
cross-sections of the TPDI of He is limited. Only time-dependent close-coupling (TDCC) 
calculations have been reported at photon energies close to double ionization threshold 
(Colgan k Pindzola 2002, Hu et al 2005). In the meantime, the fully-resolved triple- 
differential cross-section (TDCS) contains the most detailed information on the two- 
electron break-up process. In the case of one-photon double ionization, the analysis 
of the TDCS brought a wealth of information allowing to clearly separate various 
mechanisms of photo double ionization (Briggs k Schmidt 2000, Knapp et al 2002) 

So there is a compelling motivation to perform further studies of He TPDI TDCS 
at a wide photon energy range. To serve this purpose, we apply the convergent close- 
coupling (CCC) formalism to describe the two-electron continuum following the TPDI 
of He. This method has been extensively used to study one-photon double ionization 
(Kheifets k Bray 1996, Kheifets k Bray 1998a, Kheifets k Bray 19986). In the two- 
photon case, a further development of the method is needed. Firstly, integration over all 
the intermediate states of the target is required following absorption of a single photon. 
This involves evaluation of the ill-defined continuum-continuum dipole matrix elements. 
To circumvent this difficulty, we perform our calculations in the Kramers-Henneberger 
gauge of the electromagnetic field (Ivanov k Kheifets 2005 b). As an alternative and 
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much less time consuming method, we use the closure approximation to carry out 
summation over all the intermediate target states. As the result of this procedure, 
we end up with evaluation of the monopole and quadrupole matrix elements between 
the correlated ground state and the CCC final state. This procedure has been worked 
out in an earlier paper on the second Born treatment of the electron impact ionization- 
excitation and double ionization of He (Kheifets 2004). 

Secondly, as compared to the weak-field single-photon double ionization, the 
theoretical description of TPDI requires an account of the field modification of the target 
states. Conceptually, this procedure has been developed in a way which can incorporate 
the CCC formalism (Ivanov & Kheifets 2005a). However, at the present stage, we are 
unable to perform fully converged numerically accurate calculations of this type due 
to limitations of computer power. Leaving this task for the future, we report in the 
present paper the CCC calculations of the He TPDI TDCS in the perturbative regime 
when interaction of the target and the electromagnetic field is treated to the second 
order of the perturbation theory. Although this type of calculation would be clearly 
inadequate to obtain accurate total TPDI cross-sections, it can be useful in evaluating 
the angular correlation pattern in the two-electron continuum. In this respect, the 
present model is somewhat analogous to the asymptotic Coulomb treatment of the one- 
photon double ionization (Schwarzkopf et al 1994, Maulbetsch & Briggs 1993). The 
asymptotic three-body Coulomb (3C) wave function is only correct at large distances 
from the nucleus. The fact that calculations with this wave function reproduce correctly 
the angular correlation pattern in the two-electron continuum tells us that it is the 
large distances which are responsible for forming this pattern. In the meantime, the 
magnitude of the double photoionization cross-sections calculated with the 3C wave 
function is typically off by a significant factor which means that the photoelectron flux 
originates from vicinity of the nucleus where the 3C wave function is incorrect. The same 
benefits and disadvantages of such an approximate model can be seen in the present 
work. Indeed, the fact that the angular correlation pattern of the two-photon double 
ionization can be calculated correctly within a lowest order of the perturbation theory 
tells us that it is the inter-electron correlation in the final state that is responsible for 
forming this pattern. In the meantime, the total cross-section which is incorrect in such 
a perturbative calculation is sensitive to the detail of the electromagnetic interaction 
and requires a fully non-perturbative treatment. 

The paper is organized as follows. In Sec. 2 we give an outline of the formalism. In 
Sec. 3 we give our numerical results which are obtained both with the CCC integration 
over the intermediate target states and the closure approximation. Subsections 3.1, 
3.2, 3.3 contain the analysis of the integrated and differential cross-sections and the 
symmetrized ionization amplitudes, respectively. We conclude in Sec. 4 by outlining 
our future directions in this project. 
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2. Formalism 

2.1. Second order perturbation theory 

We use the following second order perturbation theory expression for the TPDI TDCS: 

> /> ^ (^/(fci)! e • d |^(p))<^(p)| e • d |^q> ' 
<J dp Eo + u-^-ei + iS {2lJ 

Here the two-photon ionization constant C 2l = 8ir 3 c~ 2 aQT = 2.505 x 10~ 52 cm 2 s" 1 (Tang 
& Bachau 1993). The vector e represents a linearly polarized light. Generalization for 
an arbitrary polarization is straightforward (Manakov et al 1996). The dipole operator 
d = d\ + d 2 where d a = r a , V Q and Zr a jr\ in the length, velocity and acceleration 
gauges, respectively, the nucleus charge Z = 2 for helium. The exponent (3 depends on 
the gauge of the electromagnetic interaction, (5 = 2 in the length gauge. 

In the CCC formalism, we represent the two-electron state by a close-coupling 
expansion over the channel states each of which is composed of a target pseudo state / 
and a Coulomb wave k: 

Here (kf\T\jk') is half-on-shell T-matrix which is found by solving a set of coupled 
Lippmann-Schwinger equations (Bray & Stelbovics 1995). The final state with two 
electrons in the continuum in Equation (1) is obtained by projecting the positive energy 
pseudostate of the matching energy Sf = k 2 /2 on the Coulomb wave \k 2 ) . 
Using the following partial wave expansions for the T-matrix 

(kf\T\jk>)= £ ^^^.^T^mH^VK) (kLnflfWTjWn^k'L') 

L,L',J 
M,M',Mj 

and the dipole matrix element: 

(kf\ e d |tf > = E e M P E^'^ ! Y lm {n) ( lf 1 1 ) (klriflfW d ||tt > 
m p im \m f M P m J 

we can transform Equation (1) into the form: 



dVtidVt 2 dE 2 



E T,{e®e}j-y™(n in2 ) (3) 

.1=0,2 hh 



X 



)h+he l[S h {k l)+ 5 l2{ k 2 )] V^jik^h) 



Here we introduced a bipolar harmonic (Varshalovich 1988) 

y l j M (nun 2 ) = {Y h { ni )®Y l2 {n 2 )} JM = ]T C7 I ^ >Iafna y, imi (n 1 )y fama (n 2 ) (4) 



m\m2 
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The unit vectors = fej/fcj are directed along the photoelectron momenta. The tensor 
{e®e} j M} = J2m p m' p CiMp im' e M P ^M' p represents the polarization of light. The reduced 
matrix element of the TPDI is given by the expression: 

u J Hiq + u — p 1 1 — Si + 10 

Here the CCC two-electron state (2), stripped of its angular dependence, is defined as 
* /l(J) (Jfe) = \\kl riflf) + E|^ f^#^ Wk'l'njlj) (6) 

jL' I 

The bare dipole matrix element between the correlated He atom ground state and the 
CCC channel state (klnili\\d\\^/o) is evaluated elsewhere (Kheifets & Bray 1998c). By 
similar technique, the bare dipole matrix element between two CCC channel states 
(kih riflf\\d\\nilipl) breaks down into one-electron radial integrals and simple angular 
coefficients. 

Expression (6) contains the dipole matrix elements between two continuum 
states. This matrix elements are ill-defined in the length and velocity gauges of the 
electromagnetic interaction. To deal with these integrals, we use the so-called Kramers- 
Henneberger form of the Hamiltonian describing interaction of the atom and the 
electromagnetic field. The matrix elements of the electromagnetic interaction operator 
in this representation can be written as (Ivanov & Kheifets 20056): 

Z Z 



a,n + p 



rVKH 
^int 




fo, n) — — — 2, / cospO 



r\ \ri + F cos9/u 2 \ 



b) d9 (7) 



We employ here the notation \a,m) where a stands for a set of quantum numbers 
describing the atom and m denotes a number of laser photons in a given mode. As can 
be seen, matrix elements defined by the Equation (7) are finite and well-defined even if 
both a and b are continuum atomic states. In Eq.(7) F is the classical field strength. 
It is a vector directed along the polarization vector of light, its magnitude is related to 
the photon density as F 2 /8tt = nuo. Both matrix elements in Equation (1) correspond 
to the dipole transitions with p — 1. For such a dipole transitions, in the weak field 
limit F — > 0, operator (7) coincides, with the dipole operator in the acceleration gauge 
Zri/rf. In addition, operator (7) can connect directly the initial and the final states in 
Equation (1) by simultaneous absorption of two photons (p = 2). Both sequential and 
simultaneous absorption processes are of the same order and should be included in the 
second order perturbation theory. See (Ivanov & Kheifets 20056) for detail. 

If we assume that the integrand in Equation (5) is a smooth function, we can take 
out an average energy denominator and use the completeness of the CCC basis. This 
procedure, known as the closure approximation, will take us to the following result: 

V^jikiM) = A- 1 (^ 1/(J) (A; 1 )|| dd \\V )(l 2 k 2 || /) , (8) 

The reduced matrix elements entering Equation (8), in the length gauge, have been 
worked out in an earlier paper on the second Born treatment of the electron impact 
ionization-excitation and double ionization of He (Kheifets 2004). 
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By integrating TDCS (3) over the angles fli, Q 2 one gets the single differential, with 
respect to the energy, cross-section (SDCS). In the closure approximation, it is given by 
the following expression: 

da ~ 2 ' -c^tAw E IcSiorEl^w)^)!!^!!^)^!!/) 2 (9) 



dE 2 7 A 2 2VE2 J=0 y 
Implicit in the above is that E\ + E 2 = Eq + 2u with E 2 = e n2 i 2 for some pseudostate 
n 2 and every l 2 . Variation of the Laguerre exponential fall-off for each l 2 allows this for 
one value of E 2 (Bray & Fursa 1995). 

Further integration of Equation (9) over the energy E 2 leads to the total integrated 
TPDI cross-sections (TICS): 

2 

= £ \cZo\ 2 E\^fL(J)(k)\\ dd ||* >| 2 (10) 

^ J=0,2 fL 

In this expression, the energy of the pseudostate and the Coulomb wave are bound by 
the energy conservation e/ + k 2 /2 = 2uo + E . When evaluating expressions (9) and 
(10), we can assume that the largest contribution to the sum over the intermediate 
states in (5) comes from those terms in which the energy of the continuum electron 
in the intermediate state is of the order of the energy of the ejected photoelectron 
p 2 /2 + ti ~ kj/2 + tf = E + 2u. Therefore we can aprroximate the ratio by A 2 /u 2 ~ 1. 

2.2. Two-photon double ionization amplitudes 

Although expression (3) has been derived in the second-order perturbation theory, 
its tensorial structure is general for the TPDI. It allows to introduce a simple 
parametrization of the TDCS in the manner suggested for the one-photon double 
ionization by Huetz and co-workers (Huetz et al 1991, Malegat, Selles & Huetz 1997, 
Malegat, Selles, Lablanquie, Mazeau & Huetz 1997). We proceed as follows. We rewrite 
Equation (3) as 



d z a 2l 

oc 



dQidQ 2 dE\ 



^6 -% W2 Kn 2 )M Zl/2 (A; 1 ,A; 2 ) (11) 



+ b 2 -y 2 L ^(n l n 2 )Q hh (k u k 2 ) 



2 









(12) 



Here we introduced the monopole M^ 2 (A;i, k 2 ) and quadrupole Qhi 2 (ki, k 2 ) reduced 
matrix elements modified by the phase factors and overlaps. The tensor bj = {e <g> e}j 
represents the polarization of light. 

We first deal with the quadrupole amplitude 



a q= E+ E + E W-y^Km^Qi^hM) (13) 

I l 1= l 2 h=o l 2 =0 ) 

= J2 b 2-y2 U (ri 1 n 2 )Q ll (k 1 ,k 2 ) 

! (n 1 n 2 )QH+2(A;i, k 2 ) + y 2 L+2 ( (n 1 n 2 )g i+2 k 2 ) 



1=0 
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We introduce the symmetrized quadrupole matrix elements: 

QthihM = \{QhiAki,k 2 )±Q hl2 {k 2l k l )} = ±Qf ih (k 2 ,h) (14) 

Using this notation and the symmetry property 3^j +2 '(n 1 n 2 ) = yj +2 (n 2 rii), we can 
write the quadrupole amplitude as 

A Q = 5> 2 - y 2 ll (n 1 n 2 )Q ll (k 1 , k 2 ) (15) 



1=0 



+ E Qhi 2 (ki, k 2 ){b 2 • y 2 " +2 (nin 2 ) + b 2 ■ y 2 ll+2 (n 2ni )} 
1=0 

+ E Qhi 2 (ki, k 2 ){b 2 ■ y 2 ll+2 { ni n 2 ) - b 2 ■ y 2 ll+2 (n 2ni )} 

1=0 

We use the bipolar harmonics expressions given by Manakov et al (1996): 

y l 2 \n u n 2 ) = Ci^P[{x){ ni <g> n 2 } 2 + P"(x){[ni x n 2 ] <g> [m x n 2 ]} 2 } (16) 

y l 2 l+2 ( ni , n 2 ) = D^Pl'ix^m ® m} 2 + P{' +2 {x){n 2 <g> n 2 } 2 (17) 

-2P;' +1 (x){n l ®n 2 } 2 } 

where x = cos 9i 2 = n\ ■ n 2 and Ci, Di are the normalization coefficients. These 
expressions allow for the following parametrization of the quadrupole amplitude: 



■b 2 



(18) 



(19) 



A Q = g + [{nj <g> ni} 2 + {n 2 ® n 2 } 2 ■ b 2 + g~ [{n x <g) n 1 } 2 - {n 2 <g> n 2 } 2 
+ &{ni <E> n 2 } 2 ■ b 2 + 5r {[^i x n 2 ] <g> [rii x n 2 ]} 2 • b 2 
Here we introduced the four symmetrized amplitudes: 

gHkuk^x) = Y.D l [P['{x) + P[' +2 {x)]Qf l+2 {k l ,k 2 ) 

1=0 

g s (h,k 2 ,x) = J2CiPl(x)Qu(k 1: k 2 )-AD l P; , +1 (x)Qt l+2 (k 1: k 2 ) 
i=i 

g (h,k 2 ,x) = ^C l P; , (x)Q ll (k 1 ,k 2 ) 

1=2 

For equal energy sharing g~{k\ = k 2 ) = 0. Equation (18) allows to separate the 
kinematic variables represented by the tensor products of vectors rii and n 2 and the 
dynamics of two-electron escape which is contained in the TPDI amplitudes. 

For the coplanar geometry, when the two photoelectron momenta and the photon 
polarization vector belong to the same plane, we can make the following choice of the 
coordinate frame: z \\ k\ and y \\ [ki x k 2 ]. In this case the quadrupole amplitude is 
given by the formula 



3* 



P 2 (cos0i) + P 2 (cos0 2 )] +g-[P 2 (cos9 1 ) - P 2 (cos9 2 )\ (20) 

1 1 ) 

+ -g s [3 cos 6 1 cos 9 2 - x\ + -g (x 2 - 1) > 
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where cos#j = rii ■ e, x = cos(#2 — The monopole amplitude is given by the 
expression: 

fo = E &o-% W Hn 1 n 2 )M W2 (fc 1 ,fc 2 ) = -4=E(- 1 ) Z ^?^ P '( x ) M »(^i^ 2 ) (21) 
h=h v3 /=0 4vr 

So, in case of an arbitrary energy sharing, the theoretical description of the TPDI TDCS 
requires five amplitudes as compared with only two amplitudes needed to describe the 
one-photon double ionization TDCS. This reflects a much richer dynamical structure of 
the two-electron continuum 

3. Numerical results 

3.1. Integrated cross-sections 

We start presenting our numerical results with the total integrated cross-section (TICS) 
of the TPDI of He. On the left panel of Figure 1 we show the TPDI TICS from several 
closure calculations at various photon energies in comparison with the literature values. 
As is seen from the figure, the presently calculated cross-section falls substantially below 
predictions of other methods. We cannot say whether this is a shortcoming of the 
closure approximation or the second-order perturbation theory in general. The full 
CCC calculation of the total cross-section is too much time consuming to get the data 
across a wide photon energy range. 




40 42 44 46 48 50 52 2 4 6 8 10 12 2 4 6 8 10 12 

Photon energy (D (eV) Electron energy E[ (eV) Electron energy Ej (eV) 

Figure 1. Left panel: Total cross-section of the TPDI on He as a function of the photon 
energy ui. Present CCC closure calculation is compared with the TD-calculation of 
Piraux et al (2003), i?-matrix calculation of Feng & van der Hart (2003) and Grid 
calculation of Hu et al (2005). Central panel: Single differential cross-section of the 
TPDI of He at to = 45 eV. The extrapolated CCC SDCS at equal energy sharing is 
shown by an open circle. The TDCC results are from Colgan & Pindzola (2002). Right 
panel: Single differential cross-section of the single-photon double ionization of He at 
the same excess energy. A "raw" and extrapolated CCC calculations are shown by the 
red solid and purple dashed lines, respectively. 

Evaluation of the magnitude of the TDCS requires the knowledge of the single- 
differential cross-section (SDCS) resolved with respect to the photoelectron energy. As 
was shown in the case of the single-photon double ionization, the "raw" SDCS, as 
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calculated by the CCC method, is oscillatory due to explicit distinguishability of the two 
photoelectrons. A special extrapolation procedure was devised to cure these unphysical 
oscillations which relied on the fact that the area beneath the SDCS is the total cross- 
section and the equal energy SDCS should be given by the coherent sum of the direct and 
exchange CCC amplitudes (Kheifets & Bray 2002). The raw TPDI SDCS from the CCC 
closure calculation at u = 45 eV is shown in the central panel of Figure 1 in comparison 
with the TDCC calculation of Hu et al (2005). The full extrapolation procedure was not 
implemented in the case of TPDI except for the midpoint which was calculated with the 
coherent sum of the direct and exchange amplitudes. For comparison, in the right panel 
of Figure 1 we show the SDCS of the single-photon double ionization of He at the same 
excess energy. Both the raw and extrapolated calculations are shown which resemble 
respectively the CCC and TDCC TPDI results on the central panel of Figure 1. 

3.2. Triple differential cross-sections 

The full CCC calculation of the TDCS TPDI based on Eqs. (3) and (5) is very time 
consuming. So far, we were able to perform just one such calculation at the photon 
energy of 45 eV and equal energy sharing of the two photoelectrons E\ = E 2 = 5.5 eV. 
The resulting TDCS are shown in Figure 2 in comparison with the TDCC calculation 
of Colgan & Pindzola (2002). Here the coplanar geometry is assumed with one electron 
escaping at a fixed angle Q\ and the second electron detected on the full angular range. 
The CCC basis of the final state included 17 — / target states with orbital momentum I 
ranging from to 4 (the so-called 17/4 calculation). The CCC basis of the intermediate 
states was much shorter with only three target states (Is, 2s, 2p) included. This was 
possible because the absorption of a single 45 eV photon could populate only few target 
states. 




90 180 270 360 90 180 270 360 90 180 270 360 90 180 270 360 

Variable angle 9t (deg) 



Figure 2. TDCS TPDI of He for the coplanar geometry at u> = 45 eV and 
Ei = E 2 = 5.5 eV. A 17/4 full CCC calculation (red solid line) and a 20/5 closure 
calculation (blue dashed line) are compared with the TDCC calculation of Colgan & 
Pindzola (2002) (black thick solid line). 

As it will be shown in the subsequent analysis, the TDCS is dominated by the 
D-w&ve except for the special case of the fixed electron angle Q\ = 90°. In this case 
a very small TDCS is a result of a delicate compensation of the D and 5-waves. The 
CCC calculation reproduces the TDCC calculation of Colgan & Pindzola (2002) fairly 
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well except for this special case when the D wave remains uncompensated and produces 
a very large peak. 

The closure approximation (8) makes a considerable simplification of the calculation 
and a much larger basis of the final target states can be handled with the present 
computational resources. In Figure 2 we show the results of a 20Z5 closure calculation 
at the same photon energy u = 45 eV and E\ — E 2 — 5.5 eV. At 9\ = 0°, both the full 
CCC calculation and the closure calculation perform fairly well, the closure calculation 
resulting in a somewhat lesser "wings" on the sides of the main peak. At 9\ = 30° 
and especially at 6*i = 60°, the full CCC calculation reproduces better a second peak 
appearing at around 62 — 270°. In the most difficult case of Q\ = 90°, both calculations 
produce a large spurious peak at 82 — 270° which is not visible in the TDCC calculation. 
As will be discussed in the following, this peak is a result of an undercompensated D 
wave. Generally, the closure approximation results are quite close to the full CCC 
calculation. In the rest of the paper, we will focus on the closure approximation which 
is much less computationally demanding. 



7 1 — . — . — 1 — . — . — 1 — . — . — 1 — . — . — 1 5 




90 180 270 360 90 180 270 360 90 180 270 360 90 180 270 360 




Variable angle 9 2 (deg) 

Figure 3. TDCS TPDI of He for the coplanar geometry at u — 42 eV and 
Ei = E2 = 2.5 eV. A 20Z6 closure calculation (red solid line) is compared with the 
TDCC calculation of Hu et al (2005) (black thick solid line). The top row shows the 
TDCS with combined contributions of D- and S*-waves whereas in the middle and 
bottom rows the separate contribution of the D and iS-waves are plotted. 

Further analysis of the TPDI TDCS is carried out in Figure 3 where we present 
a set of TDCS for the photon energy u = 42 eV and E\ = E 2 = 2.5 eV. Following 
Hu et al (2005), we give a separate contribution of the D-and S- waves to the TDCS at 
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various fixed angles 0\. As is seen in the figure, the .D-wave contribution is dominant 
for all fixed electron angles except for 9\ = 90°. Here a very small cross-section of Hu 
et al (2005) is result of an almost complete compensation of the partial S- and D-wave 
contributions. The present CCC closure calculation deviate insignificantly from Hu et al 
(2005) but no such a perfect compensation of the partial waves occurs. As the result, 
the TDCS at 6\ = 90° is quite different from the TDCC calculation. 



6 i . . i — — i . . j 6 i . . i — i — . i . . — i — 6 




Variable angle 9t (deg) 

Figure 4. TDCS TPDI of He for the coplanar geometry at w = 46 eV and various 
energy sharings. A 20/5 closure calculation (red solid line) is compared with the TDCC 
calculation of Hu et al (2005) (black thick solid line). In all panels, the fixed electron 
angle is B x = 30°. 

More data is presented in Figure 4 where we show the TDCS at the photon energy 
uj = 46 eV and various, both equal and unequal, energy sharings. Here the fixed 
electron angle is kept constant at Q\ = 30°. Generally, a good agreement with the 
TDCC calculation of Hu et al (2005) is seen with some minor shape variations. This 
has been the case previously as at this particular value of the fixed electron angle the 
main contribution to the TDCS comes from the .D-wave and no significant interference 
between the partial waves occurs. 

As it was pointed out when discussing the SDCS calculations, the equal energy 
sharing TDCS can be calculated fully ab initio in the CCC method by constructing 
the coherent sum of the direct and exchange ionization amplitudes. Therefore, all the 
TDCS shown in Figs. 2-3 are absolute. When making comparison with the TDCC 
results, only one scaling constant was used across all the panels of these figures to 
account for somewhat different total integrated cross-sections in the two methods. For 
the asymmetric energy sharing, however, the magnitude of the CCC TDCS can be 
affected by the SDCS oscillations. Therefore, the CCC results in Figure 4 were scaled 
individually on the side and central panels to make a shape comparison with the TDCC 
data of Hu et al (2005). 

3.3. TPDI amplitudes 

The TDCS presented in the preceding section can be most conveniently analyzed using 
parametrization (20) with the symmetrized amplitudes defined by Equation (19). In 
the case of equal energy sharing, only three quadrupole amplitudes and one monopole 
amplitude are needed. The moduli squared of these amplitudes are shown in Figure 5 
for uj = 45 eV and E x = E 2 = 5.5 eV. 
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Figure 5. Symmetrized amplitudes g + , g s , go and fo as functions of the mutual 
electron angle 6*12 for 10 = 45 eV and E\ = E 2 — 5.5 eV. 



As in the case of the single-photon double ionization, all the amplitudes peaks 
strongly at 812 = 180°. This is a result of the strong electron repulsion favoring back- 
to-back emission. The central part of the amplitudes around #12 = 180° was fitted with 
the Gaussian on a flat pedestal: 

\g\ 2 oc aexp 

The fitting seems to be nearly perfect with the width parameters A# 12 of 110°, 131°, 
63° and 105° for g + , g s , g and f , respectively. For comparison, at the present 
energy sharing, the symmetric amplitude of the single-photon double ionization has 
the Gaussian width of 96°, somewhat less than that of g + . 

Gaussian ansatz (22) combined with Equation (20) allows to analyze the evolution 
of the TPDI TDCS as a function of the fixed photoelectron angle. At 9\ = 0, 6u = #2 
and 

P2(cos6'i) + P2(cos6 ) 2) = -(3COS6 12 + 1) , - (3 cos 6*i cos 6*2 — x) = cos 6*2 

We see that the kinematic factors accompanying the largest amplitudes g + and g s both 
peak at 9 2 = 180° where the amplitudes have the maximum. This produces a bold 
peak seen in all TDCS figures at this combination of angles 9\ and $2- This is in sharp 
contrast to one-photon double ionization in which the kinematic factor is represented 
by Pi (cos #1) + Pi (cos #2) which has a node at #12 = 180°. As a result, the one-photon 
TDCS has a maximum at a compromise angle where neither the kinematic factor nor 
the amplitude have their respective maxima. 

The shape of the TPDI TDCS changes completely at 9i = 90° where 

p2(cos 6>i) + p2(cos 9 2 ) = -(3 cos 2 9 2 — 2) , -(3 cos^i cos #2 — x) = — - sin 9 2 

Neither of the factors peak at 180° and the P-wave contribution has the maximum at 
a compromise angle of 270°. A partial overlap between the kinematic and dynamic 
factors result in a much smaller cross-section. Coincidentally, the S"-wave has the peak 
at 9 2 = 270° as well. However, the relative phase of the S and D amplitudes conspires 
to cancel out the two contributions. This cancellation is a dynamic property depending 
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on the relative magnitude of the quadrupole and monopole amplitudes. It is somewhat 
different in the CCC and TDCC calculations which assume different models of the two- 
electron dynamics. 

In Figure 6 we analyze the TPDI amplitudes much further away from the double 
ionization threshold at the equal energy sharing of E\ = E% = 20 eV. All the amplitudes 
were fitted with the Gaussian ansatz. The width parameters are 94°, 131°, 90° and 90° 
for g + , g s , g and / respectively. There is no systematic change of the Gaussian width 
parameters at 40 eV as compared to 4 eV excess energy. In the meantime, the Gaussian 
width of the single-photon ionization amplitude has increased from 96° to 103°. 

Qualitative similarity of the amplitudes in Figure 5 and Figure 6 allows us to predict 
that the basic shape of the TDCS will remain unchanged much further away from the 
double ionization threshold compared to what was investigated by Hu et al (2005). 
To make more quantitative predictions for the TDCS a more systematic study of the 
amplitudes is needed including a thorough test of convergence with respect to the size 
of the CCC basis. 




90 180 270 360 90 180 270 360 90 180 270 360 90 180 270 360 

Mutual angle 6 12 , deg 

Figure 6. Symmetrized amplitudes g + , g s , go and fo as functions of the mutual 
electron angle 612 for E\ = E 2 = 20 eV. 

4. Conclusion 

In the present paper we investigated the two-photon double ionization of He in a 
wide range of photon energies from 4 to 40 eV above the threshold. The electron- 
photon interaction was treated to the lowest second order whereas the electron-electron 
interaction was included non-perturbatively. We applied the convergent close-coupling 
formalism to describe the S- and D two-electron continua of the final doubly ionized 
state. The intermediate P-state was also described with the CCC expansion, albeit of 
a smaller size. The ill-defined continuum-continuum dipole matrix elements between 
the intermediate and final states were handled using the Kramers- Henneberger gauge 
of the electromagnetic interaction. As an alternative and much less computationally 
extensive solution, we used the closure approximation summing over the complete set 
of the intermediate states. By doing so, we were able to reduce the TPDI amplitude 
to the matrix element of the squared dipole operator, in the length gauge, between the 
correlated ground state and the two-electron CCC final state. 

Near the threshold, the presently calculated total integrated TPDI cross-section is of 
the same order of magnitude as predicted by other non-perturbative methods. However, 
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further away from the threshold, the CCC closure results fall substantially below the 
reported literature values. The single differential, with respect to the energy, cross- 
section was also evaluated. As in the single-photon double ionization, the raw CCC 
SDCS shows unphysical oscillations and needs rescaling. However, the equal energy 
sharing point can still be calculated ab initio. 

The fully-resolved triple differential cross-sections were evaluated using integration 
over the complete set of CCC intermediate state as well as the closure approximation. 
Both sets of results agree reasonably well with the previously reported TDCC cross- 
sections. The only exception is the fixed electron angle of 90° where the delicate 
cancellation of the S- and D-partial waves takes place. Here the CCC and TDCSS 
results are at variance. These findings indicate that the angular correlation pattern in 
TPDI is formed mainly due to the inter-electron interaction in the two-electron continua. 
Influence of the strong laser field and precise mechanisms of the two-photon absorption is 
less important. Of course, this conclusion needs to be further confirmed experimentally. 

A complete set of symmetrized amplitudes was introduced to parametrize 
conveniently the TPDI TDCS. The amplitudes display a Gaussian shape and allow 
to explain the evolution of the TDCS at varying fixed electron angles. The amplitudes 
evaluated at fairly high excess energy of 40 eV resemble closely the amplitudes near 
threshold. This allows to predict a rather robust shape of the TPDI TDCS across a 
wide range of photon energies. Reported Gaussian parameters can be useful in modeling 
the TPDI TDCS at various experimental kinematics. A more systematic and consistent 
study is required, and is planned, to serve this purpose. Numerical values reported 
here serve only as an illustration and should be treated as preliminary because not all 
the convergence checks were possible to perform at the present stage. Therefore, some 
minor variation in the width parameters might occur as the size of the CCC basis varies. 

This work is the first report on the application of the CCC method to the TPDI 
process in He. We intend to continue this work, both within the scope of the perturbation 
theory and treating the electromagnetic field non-perturbatively. The basic framework 
of the theory is outlined in our earlier paper (Ivanov & Kheifets 2005a). The theory 
is potentially capable of studying a wide range of photon energies and electromagnetic 
field intensities. The only limitation is available computational resources which we hope 
to lift in the near future. 
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